Out-of-equilibrium dynamical fluctuations in glassy systems 
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In this paper we extend the earlier treatment of out-of-equilibrium mesoscopic fluctuations in 
glassy systems in several significant ways. First, via extensive simulations, we demonstrate that 
models of glassy behavior without quenched disorder display scalings of the probability of local 
two-time correlators that are qualitatively similar to that of models with short-ranged quenched 
interactions. The key ingredient for such scaling properties is shown to be the development of 
a critical-like dynamical correlation length, and not other microscopic details. This robust data 
collapse may be described in terms of a time-evolving Gumbel-like distribution. We develop a 
theory to describe both the form and evolution of these distributions based on a effective sigma- 
model approach. 



PACS numbers: 61.20.Lc, 75.10.Nr 
I. INTRODUCTION 



The dynamics of glassy systems is extremely slow, 
both near and far from equilibrium 1 . Until recently, ex- 
periments, numerical studies and theory focused on the 
study of averaged global, or bulk, dynamical quantities 
that demonstrate the nonequilibrium character of the 
glassy state. Interest is now shifting to the study of local 
quantities embodied in intermittent fluctuations and dy- 
namic heterogeneities. The recent development of pow- 
erful experimental techniquesS^^SiLSiSiiSi and the use 
of extensive numerical simulations have started to yield 
valuable information about the dynamics of supercooled 
Uquida 11 ! 12 ! 13 ! 14 ! 15 . 16 ! 17 ! 18 ! 19 and glassy systems^! at a 
mesoscopic scale. A detailed analysis of temporal fluctu- 
ations in the evolution of global observables allows one to 
distinguish between intermittent and continuous dynam- 
ics, the former being associated with sudden, large re- 
arrangements&Siifi. In addition, the study of run-to-run 
fluctuations in systems of mesoscopic size yields comple- 
mentary information 21 . The insight obtained from these 
studies is useful in identifying the most relevant rear- 



rangements that take place during evolution and thus 
may lead to a better understanding of the mechanism 
leading to glassy arrest. 

The behavior of global dynamic correlations in macro- 
scopic glassy systems is rather well-described analytically 
with an extension of the mode-coupling theory to the low 
temperature regim o 22 : 23 . This approach is also attrac- 
tive since, via its connection to mean-field models with 
quenched disorder, it allows for an interpretation of dy- 
namic arrest and slow dynamics in terms of the geomet- 
ric properties of a free-energy landscape 2 ^ 2 ^ 2 ^!. How- 
ever, it is not clear how to describe temporal and spatial 
fluctuations within this approach. Another picture of 
glassy dynamics is based on kinetically constrained mod- 
els with trivial static interactions2Si22^£ and is successful 
in describing aging effects and other glassy features^ as 
has been recently reviewed 3 ^. Moreover, dynamic hetero- 
geneities also find a natural interpretation within these 
model o 33 ' 34 . The connection between the microscopic in- 
teracting real system and the non-interacting kinetically 
constrained model is however phenomenological, as the 
facilitation rules are based on a coarse-grained descrip- 
tion of the liquid dynamics. Other phenomenological de- 
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scriptions do not even have a real-space interpretation. 
Thus, one of the more interesting and open theoretical 
problems at present is how to derive the observed proper- 
ties of medium-size and large-scale fluctuations that are 
possibly associated with cooperative and correlated re- 
arrangements of the microscopic constituents (particles, 
spins, etc.) from a microscopic model 3 ^. 

In this paper we focus on particular aspects of the non- 
equilibrium dynamics near or below the onset of glassy 
arrest. Some of our arguments and numerical analysis 
can be simplified to capture the (stationary) slow dy- 
namics of supercooled liquids as a particular case. The 
purpose of this paper is three-fold. First, we test whether 
the scenario proposed in Refsi22*2ii2£ for spatially het- 
erogeneous glassy dynamics applies to models without 
quenched disorder such as kinetically constrained spin 
models. Second, we show that several features of the 
probability distributions functions (pdf's) of fluctuating 
two-time correlations, namely time scalings and func- 
tional forms, are common to many glassy systems. Third, 
we discuss an analytic approach developed to describe 
the scaling form and time evolution of these pdf's. This 
approach is based on the arguments put forward i n 20 i 35 
and on an analogy with the study of the order-parameter 
fluctuations in critical phenomena™, where one can de- 
scribe universality classes in terms of the pdf's of global 
fluctuations. The latter is similar in spirit to the work of 
Racz and collaborators on interfacial growth problems 37 , 
and Bramwell and coworkers on the xy model™. 

Even if, strictly speaking, there is no liquid-glass dy- 
namic phase transition, glassy dynamics present some 
aspects of critical dynamics. Among others, a dynamic 
correlation length defined from the fluctuations of the lo- 
cal two-time correlations shows an apparent divergence 
at long times2L2M2ilii^. In the critical situation, one 
expects that the symmetries of the problem constrain 
the form and scaling of the pdf's of fluctuating dynam- 
ical quantities. A natural way to attempt an analytic 
approach to these fluctuations is then to construct an ef- 
fective sigma model that respects these symmetries. This 
approach has been advocated inSL^LS. Here we extend 
it to propose a model that describes not only the time- 
scaling of the pdf's but also their functional form. 

The plan of the paper is as follows. In Section [Q] 
we summarize the main properties of glassy dynamics 
as displayed by the decay of global correlation func- 
tions. In Section II I II we present a short review of ex- 
periments that focus on the measurement of dynamic 
fluctuations. In Section IIVI we introduce the models 
that we study numerically; these are a modification of 
the Fredrickson-Andersen (fa) model due to Graham et 
alm& and the Edwards-Anderson (ea) spin-glass™ In 
Section we present the pdf's for fluctuating two-time 
correlations. We demonstrate that very similar pdf's are 
obtained from facilitated Ising models without quenched 
random interactions and the disordered spin-glass prob- 
lem. These forms also resemble strongly the experimental 
data of Cipelletti et al~ for colloidal suspensions. We 



describe these pdf's with an evolving Gumbel-like form 
that we later argue can be justified from the sigma- model 
approach. In Section IVII we introduce the theoretical 
modeling of these two-time dependent pdf's. Finally, in 
Section fVIII we present our conclusions. Appendix [S] dis- 
cusses the definition of a Gumbel distribution and some 
of its properties. The aim of Appendix [5] is to identify 
which features of the PDFs are trully due to glassy dy- 
namics. To this end we analyse the stochastic dynamics 
(with no facilitation rule) of a finite system of indepen- 
dent spins on a lattice. 

II. NONEQUILIBRIUM RELAXATION IN 
GLASSY SYSTEMS 

Glassy systems show slow nonexponential dynamical 
relaxation. This is demonstrated by the extremely slug- 
gish decay of one-time observables like the energy density, 
the structure factor, etc. In many experimental cases, 
one-time observables can be considered to be approxi- 
mately constant within the time window explored. Thus, 
observables that depend on only one time are not suffi- 
cient to characterize the non-equilibrium relaxation of a 
glassy system in its full richness; they do not signal one 
of the hallmarks of glassy systems, namely aging phe- 
nomena. Aging is partially described by the decay of 
two-time global correlation functions: 

C(t,t w ) = N- l {0{t)0(t w )) (1) 

with 0(t) the global observable of interest and a two- 
time dependent normalization that enforces C(t,t) — 
1. The angular brackets denote an average over differ- 
ent thermal histories of the system. In a model with 
quenched disorder, e.g. the EA spin-glass, one averages 
the global correlation over different realizations of the in- 
teractions. Sometimes a coarse-graining in time using a 
short time window around the observation instant is also 
implemented. It is important to note that for a suffi- 
ciently large system, the quantity 

0{t)0(t w )/N (2) 

does not fluctuate and the result obtained for a single 
thermal history (and a single realization of disorder) co- 
incides with the averaged one, as long as times are long, 
but not so long as to exhibit recovery of equilibrium be- 
havior. 

As already mentioned, one-time observables evolve 
very slowly in the glassy phase of real systems. Still, 
some models used to capture glassy dynamics have a 
faster evolution of one-time quantities. In these cases, 
it is convenient to study the evolution of the connected 
correlations defined as in Eq. Q where 

0(t) -> 0(t) - (0(t)) . (3) 

By the very definition of aging, older systems relax 
in a slower manner than younger ones. One defines the 
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age of a system as the time spent in the phase under 
study. A way to characterize the aging properties is to 
monitor the time evolution of the two-time correlation 
defined in Eq. Q or the two-time function in Eq. J2J. 
Well below T g (conventionally defined as the value at 
which the relaxation time reaches 100s) the system is in a 
"phase" that we call "glassy" . At such low temperatures 
the correlations depend explicitly on t w and, within the 
experimentally accessible time-window, do not show a 
tendency towards decay to equilibrium behavior, that is 
to say an independence from t w . Close to glassy arrest 
one observes "interrupted aging," namely a dependence 
on the age of the system until observations reach the 
equilibration time (t w > t EQ ) where the dynamics crosses 
over to an equilibrium one. 

The correlations of Eq. (JTJ have been studied exper- 
imentally in a number of glassy materials. In spin sys- 
tems the natural choice for 0(t) is the fluctuation of the 
global magnetization^, while in colloidal suspensions it 
is the intensity of scattered light 10 ' 46 ' 47 ' 48 . Many numer- 
ical studies of a large variety of models have focused on 
the analysis of the closely related intermediate scattering 
function. On the analytical side, the study of mean-field 
disordered models such as the p-spin model yielded sev- 
eral features of the aging behavior seen in simulations and 
experimentsSiSi. Connected global correlations in sev- 
eral kinetically constrained spin systems have also been 
studied yielding similar results^i* 3 ^. 

In structural glasses, a pictorial explanation of aging 
may be put forward by imagining that each particle sees a 
cage made of its neighbors. When r is short compared to 
a characteristic time To(t w ), each particle rattles within 
its cage, decorrelation is characterized only by thermal 
fluctuations, and correlations decay in a stationary man- 
ner from unity at equal times to a plateau, g EA (T), seen 
on a logarithmic scale. When r increases, the motion 
of the particles destroys the cages and the structure of 
the system relaxes. The waiting-time dependence implies 
that the cages become stiffer as time evolves. The motion 
of a tagged particle observed with confocal microscopy 
demonstrated this scenario 6 . Upon examination of the 
value of the displacement reached at the maximum sepa- 
ration time, however, one notices that motion occurs only 
on length scales of the order of the radius of a particle. 
This means that the structural relaxation occurs without 
large scale structural rearrangement and with essentially 
no exchange of particles within the bulk of the sample. 
The same observation applies to the dynamics of granular 
matter under external forcing^. 

The existence of a well-defined cage translates into the 
appearance of a plateau in the decay of the two-time 
correlations. The presence of a plateau in the correla- 
tion function means that there is a macroscopic time 
interval in which the correlation (when plotted in log- 
log scale) can be approximated by a (non-zero) constant. 
In this case the two-time correlation can be scaled in a 
satisfactory way by using a separation in fast and slow 



regimea 23 i 24 i 50 : 

C(t,t w ) = Cp(t — t w ) + Cg(t, t w ) 

with h F (t) = e'/ TF and r F a characteristic time scale. 
The first term describes the stationary approach to the 
plateau, the second term the ^-dependent departure 
from it. The light-scattering study of a colloidal system 
of Viasnoff et al^ shows a very clear plateau separat- 
ing the two regimes. However, a clear-cut plateau is not 
always visible in numerical simulations and experiments. 

In most cases studied so far the second (aging) decay 
is characterized by a similar form 

with h s (t) a system-dependent monotonic function. A 
typical and well-understood example in which the slow 
part of the decay is represented by this form is the case 
of systems undergoing domain growth. In this case, for 
long and well-separated times, h s (t) is the typical size of 
the domains at time t. Other examples include the relax- 
ation of the incoherent scattering function in Lennard- 
Jones mixtures and the result of light-scattering mea- 
surements in colloidal suspensionsifli^L 4 ^. Monte Carlo 
data for the two-time spin-spin correlations in the 3<i EA 
model can be described with such a two-step relaxation 
and h s (t) = Similarly, the magnetization correla- 

tions in the insulating spin-glass studied by Herisson and 
Ocio^ are described by a two-step relaxation. In this 
study, an enhanced power law, h s (t) = e ln t ^ t " , yields 
the best fit of the experimental data. We shall see below 
that, for the times that one can reach with a numerical 
simulation, the kinetically constrained models defined in 
Sect. HVl either do not show a plateau in their connected 
correlations, or show several plateaus. Moreover, it is 
very difficult to identify a single function h s (t) that de- 
scribes the slow decay. This might be due to the fact that 
the slow decay involves several correlation scales that are 
still not well separated during numerical times. 

III. HETEROGENEOUS DYNAMICS 

The existence of dynamic heterogeneities in super- 
cooled liquids and glasses has been suggested on the basis 
of experiments performed using many techniques^; two 
of them are the confocal microscopy and light scattering. 
We discuss these results in some detail since we shall 
compare our numerical simulations to them. 

Confocal microscopy allows one to reconstruct the tra- 
jectory of each particle in a colloidal suspensions made of 
several thousand particles 3 ^. Kegel and van BlaaderenA 
found a non- Gaussian distribution of particle displace- 
ments in a dense system of colloidal hard spheres and 
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they associated this large distribution to the presence of 
dynamical heterogeneities. Later, Courtland and Weeks 
found that mobile (and immobile) particles cluster in re- 
gions of lower (higher) density than average in the glassy 
phase of a colloidal system^. 

The "time-resolved light-scattering technique"^ has 
been developed with the aim of testing the fluctuations 
in time of the intensity-intensity correlations in samples 
of medium size. The existence of intermittency in mea- 
surements of global quantities may be another manifesta- 
tion of heterogeneous dynamics. Temporal fluctuations 
in colloidal suspensions have been studied in this way. 
Interestingly enough, Cipelletti's group showed that ag- 
ing samples, such as concentrated colloidal gels, have an 
intermittent dynamics leading to negatively skewed dis- 
tributions of the two-time intensity-intensity correlations 
when the time lag is shorter than the relaxation time and 
Gaussian distributions when the time lag approaches the 
relaxation time. They also showed that simpler systems 
with no glassy features, such as Brownian colloidal sam- 
ples, have a continuous relaxation with a Gaussian distri- 
bution of the same quantities for all time lags explored. 
The same trend of time resolved PDF's of particle dis- 
placements was demonstrated by Courtland and Weeks 
using confocal microscopy^. It is interesting to note that 
these experiments show that significant fluctuations oc- 
cur even if the time lag r is very short compared to t w . 

Evidence for intermittency in the global voltage signal 
noise in a solid-like colloidal glass and a polymer glass 
has also been presented recently 9 . While the samples 
are out of equilibrium the voltage time series has bursts 
with very large amplitudes. These bursts become rarer 
with the sample age in the sense that their amplitude 
decreases and the time elapsed between two consecutive 
bursts increases. When aging stops the bursts disappear 
from the noise signal. The waiting-time dependent PDF, 
p(V), constructed by dividing time in intervals of dura- 
tion r is clearly negatively skewed at short times and 
progressively approaches a Gaussian distribution when 
the age of the system increases** 

On the numerical side, most studies have fo- 
cused on the supercooled liquid phase of binary mix- 
ture9ii*i2iiii^iiSii& and polymer meltsi£ and the hetero- 
geneities have been studied by tagging each particle and 
classifying them according to their mobility during a se- 
lected time interval. Some recent papers treat the glassy 
phase of the Lennard-Jones system in a similar way^. 
Following the dynamics of individual particles is very use- 
ful to get an intuitive understanding of the evolution of 
the system, but it is hard to use as a direct input in a 
theoretical approach. In order to develop such a descrip- 
tion, it is much more convenient to use the coarse-grained 
two-time quantifie d 20 ' 35 that we shall define in Sect. [V] 
In this setting, a heterogeneity is a region in the sam- 
ple that relaxes very differently from the bulk. Clearly, 
the heterogeneities are not static and they appear and 
disappear as dynamic fluctuations. 



IV. THE MODELS 

We shall focus on several spin models that describe and 
capture different aspects of glassy systems. We define 
them below. 



A. Kinetically constrained model 

1. A variation on the Fredrickson- Andersen model 

The first model we study is a kinetically constrained 
spin system with trivial equilibrium properties. It is 
a simple variation of the Fredrickson- Andersen (fa) 
models due to Graham et al. 43 . N Ising spins, s, = ±1, 
i = 1, . . . , N, occupy the vertices of a regular cubic lattice 
in d dimensions. The energy is given by 

h N 

£ = -fE s - (6) 

Thus there are no static interactions between the spins 
but an external field of strength h is applied. The sys- 
tem is in contact with a heat bath at temperature T 
that, for convenience, we shall measure in units of the 
external field h. This model does not have a static tran- 
sition and at equilibrium the averaged magnetization, 

to = J2i( s i)/ N > is g iven b y 



The non-trivial evolution of this non-interacting system 
arises due to the constrained dynamics imposed on the 
spins. At each time step a spin, Sk, that is chosen at 
random among the ensemble is allowed to flip with a 
generalized heat bath rule. If at least p of its z = 2d 
neighbors point down, i.e. opposite to the preferred di- 
rection, the move is accepted with unit probability if 
AE < and with probability e - AE / k sT jf > 0. 
AE = Et — Ei, with Ef the total energy after the move 
and Ei the total energy before the move. In this case, 
AE = (— s£ + s|)/2 = si. Otherwise the attempted flip 
is rejected. The parameter p defines the facilitated rule. 
For instance, a three-site rule requires a minimum of 3 
neighbors to oppose the field direction. As usual, time 
units are defined via the number of attempts to flip ran- 
domly chosen spins in the sample. The facilitation rule 
satisfies detailed balance. Thus the system evolves to- 
wards the equilibrium configuration, although it may do 
so in a very slow manner. The larger the value of p the 
slower the resulting dynamics. 

Typically, one is interested in following the relaxational 
dynamics of a random initial condition in which each spin 
points in the up or down direction with probability one 
half. This mimics a rapid quench from infinite tempera- 
ture to the working temperature T. The dynamics tend 
to order the system in the up direction. However, as 
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time passes, it is harder to choose a spin with sufficient 
number of neighbors pointing down. This induces a dy- 
namic slowing down, and below T ~ 0.5 (in d = 3 and 
using p — 3) or T ~ 0.7 (in rf = 2 and using p = 2) the 
system does not manage to get close to the equilibrium 
configuration. The out-of-equilibrium relaxation below 
this crossover temperature resembles that of supercooled 
liquids^. 

Importantly, these models do not contain quenched 
disorder. The absence of quenched disorder and the local, 
non-mean-field nature of the relaxation strongly suggests 
that such models are relevant for understanding dynam- 
ics of supercooled liquids. The existence or absence of 
a genuine glass transition is not important for our pur- 
poses. Indeed, as it has been recently shown rigorously 
is some cases these models do not have a dynamic phase 
transition in the thermodynamic limilsS. Nevertheless, 
the evolution of a large finite system at long but finite 
times shares some keynote properties with many glassy 
systems, including systems that have a finite tempera- 
ture dynamic transition, such as spin-glass models that 
undergo one-step replica symmetry breaking. Thus, the 
study of facilitated models in the context presented here 
is of even broader importance than just the modeling of 
supercooled liquids. 



B. The Edwards- Anderson spin glass 

A disordered model that we shall briefly discuss here is 
the standard non-mean-field model for spin-glasses, the 
finite dimensional Edwards- Anderson model^. The en- 
ergy function is highly non-trivial in this case, 



Ej = -Y t A 

(ij) 



■8 j S ; 



(8) 



where the spins are Ising variables that occupy the ver- 
tices of a cubic d-dimensional lattice, the sum runs over 
nearest- neighbors, and the coupling strengths are ran- 
dom and chosen from a bimodal or Gaussian probability 
distribution with zero mean and variance [Jfj] = l/(2z). 
The square brackets denote an average over the random 
exchanges. This model is said to be disordered due to 
the quenched randomness introduced via the exchanges. 
The static properties are still under debate. It is be- 
lieved that in 3d the Edwards- Anderson model under- 
goes a static phase transition from a paramagnetic to a 
spin-glass phase at a finite value of T, which depends 
on the type of disorder used, but there is no consensus 
about the nature and properties of the static spin-glass 
phase'"' 4 . Dynamically, one may also follow the evolution 
of spin configurations from a random initial condition. 
The dynamics is, in this case, usual Monte Carlo with 
the heat bath rule determined by the energy function 
(JSJ. The evolution of a sufficiently large system at low 
temperatures occurs out of equilibrium. 

The local mesoscopic fluctuations of this model have 
been studied analytically and numerically nj22i2iiS. We 



shall not reproduce these results here. In this manuscript 
we merely augment the results presented in22i2iiS to test 
the properties of fluctuations that have not been dis- 
cussed in previous publications. 



V. FLUCTUATIONS 

In this section we introduce the quantities we shall 
focus on and we discuss the behavior of their pdf's as 
found for the kinetically constrained models and the 3d 
ea spin-glass. 



A. Spatial fluctuations 

Local fluctuations in an unperturbed system are mon- 
itored by studying 



C°{t 7 t w ) = i^O(i,()0(i,t,) 



(9) 



X£V r 



xGV r 



x'ev r 



where V r = £ d is a coarse-graining volume centered at 
the position r. This is a local connected function that 
evolves from a time-dependent value at equal times to 
zero at diverging time separations. Since it is convenient 
to work with a normalized two-time quantity at equal 
times, we define 



Cr (£, tw ) — 



C r (t, ty 



y/C° (t , t) C° (t w , t w ) 



(10) 



With this definition the normalization itself is site- 
dependent and might introduce spurious fluctuations. 
Another possibility is to normalize the local correlation 
with the global connected one 



C T {t 1 tw) 



C r (t : t w ) 



\J C(t, t)C(t W , ty 



(11) 



This normalization bounds the global correlation to be 
less than unity but it does not bound the local one; in 
particular, it does not enforce the local correlation to be 
one at equal times. In what follows, we present results 
obtained using definition II II only. 
For lattice spin models, the choice is 



C°(t,t w ) — — ^ K Si(t)si(t w ) 

Vr i w . 



(12) 



iev r 



V r 



jev r 



and its normalized counterpart. For the EA model one 
easily verifies that, when the coarse-graining volume 
is sufficiently large, the local magnetization vanishes, 
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K ~ 1 J2iev s i(tw) = 0, and there is no need to use the 
connected functions above. For the modified FA model, 
the local magnetization does not vanish and we shall work 
with connected functions. In none of these expressions 
is there any averaging, apart from the coarse-graining in 
space. 

With the local observables one can define a two-time 
dependent correlation length £(t,t w ) via 21 

A(A;t,t w ) = 
— / d d x 0(x, t)0(x, t w )0(x + A, t)0(x + A, t w ) 



and the normalized quantity 

A(A;t,t w ) 



B(A,t,t w ) 



A(A -> oo;t,t w ) 



A(A = 0; t, t w ) - A(A -> oo; t, t w ) 
„-A/£(t,t„) 



of the system we used only one run (there is no statistic 
average) . 



0.6 
0.5 
■ 0.4 
, 0.3 

cr 

E 0.2 
0.1 





1 



3 4 5 6 

log, (t) 



(We have made explicit the vectorial character of the 
positions £ in d > 1. In the case of spin problems it 
is convenient to use 0{x,t) = Si(t) where the overline 
indicates a coarse-graining in time over a time-window 
centered on t.) 

Depending on the relative value between the linear size 
of the coarse-graining volume, £, and the two-time depen- 
dent correlation length, £, we expect either: 

i. If £ is not much larger than £(i, t w ) different cells in 
the system behave rather differently and the histogram 
of the quantities IjlUII leads to a non-trivial PDF whose 
form depends on the two times chosen, t w and t. 

ii. If I ^> £,(t,t w ) the coarse-graining averages over 
many independent cells and the PDF approaches a Gaus- 
sian form that, in the limit £ ~ L — > oo, has a vanishing 
variance and its average is given by the global value of 
the correlation. 

As a consequence one expects to find similar fluctua- 
tions by taking a system with finite size and examining 
the behavior of the run-to-run fluctuations of the global 
two-time quantity J5J. Indeed, if L ^> £ one falls in case 
ii. Instead, if L is of the order of £, case i is found and 
the "mesoscopic fluctuations" should be non-trivial. This 
argument justifies the comparison of our approach to the 
ones in*i ancUS. 



B. Variant of the Fredrickson- Andersen model 

In this section we present the analysis of the dynamics 
of the kinetically facilitated model in d = 3 choosing 
p = 3 and then in d = 2 for p — 2. All the 3d numerical 
data have been obtained using a cubic lattice of linear 
size L = 128 lattice spacings at temperature T = 0.4 
and the 2d data using a square lattice with L = 1000 
at temperature T — 0.6. Coarse-grained data have been 
obtained using a cubic coarse-graining volume of linear 
size I = 9 lattice spacings for d = 3 and a square volume 
with £ = 20 for d = 2. For d — 3, given the very large size 



FIG. 1: Time evolution of the magnetization density in the 
modified FA model, d = 3, L = 128, T = 0.4. The line is a fit 
to the curve given in Eq. i1l:-H . 

Let us first discuss the main qualitative features of 
the evolution of the spin configuration before entering 
the quantitative analysis of the time-dependent observ- 
ables. Initial configurations are random, in that each 
spin points up or down with probability one half, which 
mimics an infinitely rapid quench from T = oo. Initially, 
approximately 66% for d — 3 (56% for d — 2) sites sat- 
isfy the facilitation constraint and are able to flip with 
the heat bath probability weight. At the early stages 
of the evolution it remains relatively easy to find spins 
that satisfy the facilitation constraint. As time evolves 
the number of positive spins, and hence the magnetiza- 
tion, increases. However, as time passes one expects to 
have isolated spins and small islands of spins that point 
downwards. These may be remnants of the initial config- 
uration or they may be due to "equilibrium-like" thermal 
fluctuations. The former are blocked until other down- 
pointing spins diffuse. The interior of the latter can flip 
but their border might be blocked due to the facilitated 
rule. A crossover (in the total time) between the "easy" 
initial regime and the "hard" latter regime was predicted 
inH. 

The analysis of the evolution for d = 3 towards its 
equilibrium value gives a first indication of the apparent 
"glass transition" temperature. Figure ^ shows the time 
dependence of the magnetization density at T — 0.4 on 
a log-linear scale. We see the above mentioned crossover 
occurring at t co ~ 10 3 MCS. The value of the magne- 
tization at the longest time reached in the simulation, 
t = 10 8 MCS, is still 10% away from the equilibrium value, 
m EQ (T = 0.4) ~ 0.848, showing that the system is still 
far from equilibrium. Assuming that the curve found for 
this range of times does indeed approach the equilibrium 
value and no further crossover to a different dynamic law 
occurs (a hypothesis that might not hold, see^ for a more 
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detailed analysis) we find that the second slower regime 
is well-characterized by the law 

= 0.4)-Alogr a i , (13) 
0.848, A = 1.37, a = 1.22 , 



m(i) ~ 

<-eq(T : 



0.4) z 



demonstrating very slow dynamics even in the behav- 
ior of this one-time quantity. At higher temperatures 
one finds that the maximum magnetization reached in 
the numerical time window gets closer to the equilibrium 
value and, for example at T ~ 0.7 it actually reaches 
this value in t ~ 10 7 time steps^. This value is consid- 
erably lower than the putative glass transition identified 
by Graham et a/. 43 . These authors might have overesti- 
mated this temperature by simulating very small systems 
(of maximum size V = 16 3 ) as compared to the ones we 
use (V = 128 3 ). 




2 4 

log, (t-t„) 

FIG. 2: The global two-time normalized connected correlation 
function, C(t,t w ), as a function of the logarithm of the time 
difference t — t w , after a sudden quench to T — 0.4. The 
curves correspond, from left to right, to seven logarithmically 
spaced values of the waiting-time, t w = 10 1 , 10 2 , 10 3 , 10 4 , 
10 5 , 10 6 , 10 7 time steps, d = 3, L = 128. 

What can we say about the structure of the remain- 
ing negative spins in the sample at long times? A set of 
snapshots of the instantaneous configuration shows that 
at around the crossover time there are no finite size do- 
mains left but only isolated spins pointing downward or, 
at most, very short strings of negative spins. It is then 
clear that their annihilation needs the cooperative dif- 
fusion of other spins in their neighbourhood, which is a 
very slow process indeed. 

The aging curves for the global connected and normal- 
ized correlation function defined in Eq. (|13fl are shown in 
Fig. |3 for several values of t w and as a function of t — t w 
in log-linear scale. As discussed in Sect. [HI disconnected 
correlation functions in glassy models usually show a two- 
step relaxation with a first rapid decay towards a plateau, 



that one estimates to be at q E A(T) 



EQ 



(T) for a model 



undergoing simple domain growth, and a second slower 
aging decay towards zero. Connected correlation func- 
tions do not necessarily show a plateau. This is the case 



in domain growth. For the times explored numerically in 
the spin-facilitated model we do not see the development 
of the plateau (which is not even seen on a log- log scale) . 
This might be due to the fact that the system is still very 
far from an asymptotic dynamic regime or it may simply 
mean that the connected correlation does not have a fast 
regime^. 
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FIG. 3: Attempt to scale the global connected and normalized 
correlation function presented in Fig. |5] See the text for a 
discussion. 

We have attempted to scale the correlation data for 
all time differences and waiting times that are equal or 
longer than t w — 10 3 time steps (note that the curves for 
the two shorter waiting times are clearly different from 
the rest). We found that the scaling 

f\n(t - t w y 



C{t,t w )~f 



\ hi t u 



(14) 



describes the data quite precisely for intermediate wait- 
ing times (t w = 10 2 -10 4 time steps), but fails to collapse 
the data for longer waiting times, as shown in Fig. [3] 
Interestingly, this scaling form is a precursor of dynamic 
ultrametricity£2i££. We have searched for other simple 
scaling forms without finding any that describe the full 
data in a satisfactory way. For this reason, we shall use 
the value of the global correlation itself to scale the PDF's 
of the local fluctuations. The same difficulty arises in the 
other kinetically constrained models that we have stud- 
ied. Henceforth we omit the subindex S that indicates 
the slow regime, since the kinetic facilitated models we 
have analyzed do not seem to have a first fast regime in 
their connected correlations. 

Let us now turn to the study of the fluctuations of 
the local correlation functions, C°(t,t w ), as defined in 
Eq. (|10f> and normalized to unity at equal times. What 
we clearly observe is that, even in this preasymptotic 
regime (in the sense that the magnetization, a one-time 
quantity, is far from its asymptotic value), the scaling of 
the PDF of local correlations is the same as the scaling 
of the global function. Moreover, the qualitative behav- 
ior of the form of the pdf's, their skewness and other 
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correlation can be described with h s (t) 
checked if 



-2 2 
(C r -<Cr>)/a 



FIG. 4: Normalized PDF of the local correlations C r (t,t w ) in 
the modified FA model, d = 3, L = 100, £ = 5 and T = 0.4. 
All curves drawn with symbols correspond to the waiting time 
t w = 10 4 MCS at various the time differences t — t w . The solid 
line represents the Gaussian normal form. The dashed lines 
are fits to modified Gumbel forms with parameter a given by 
a — 13 for t — tw = 1, a = 50 for t = 1600, and a = —13 for 
t = 400000. 



momenta, are similar to the ones that Cipelletti et al. 
found experimentally 10 . Working at different tempera- 
tures or using different linear sizes does not modify this 
conclusion. 

The summary of the features shown by the pdf's of lo- 
cal connected functions is the following. Figure 01 shows 
a first indication of the evolution of the PDF of local con- 
nected functions l|13|) for fixed t w = 10 4 time steps and 
several values of t — t w . Initially, the curves are asym- 
metric with respect to their average at short values of 
t — t w , as exhibited by the curves for t — t w = 1 and 
10 time steps. When t — t w increases they slowly ap- 
proach a Gaussian normal form as shown by the curve 
for t : — t w = 1584 time steps. Still later they go beyond 
the Gaussian form, and they skew in the opposite direc- 
tion, e.g. t — t w — 63095 time steps. Finally, the curves 
very slowly tend to approach the Gaussian form again, 
e.g. t — t w = 398107 time steps. 

The skewed pdf's resemble generalised Gumbel's dis- 
tribution (see App. IA"|) . Even if this result is most proba- 
bly just an approximation, and we cannot give it a prob- 
abilistic interpretation, it is worth mentioning the resem- 
blance here. In Sect. IVll we argue that these forms, and 
their evolution in time, can be understood using an ef- 
fective random manifold theory that naturally gives rise 
to pdf's with a Gumbel-like form. 

The next point we investigate is whether the pdf's of 
local correlations in the slow regime of widely separated 
times scale as the global correlation itself, namely 

P(C r ;t,t w ) =g(C r ;C(t,t w )) . (15) 

In22*2I we studied the 3d EA model whose slow global 



P\C r ; t, t w ) 



g(C r ; C s (t,t w )) 
'h s (t w ) 



Cr'J 



h s (t) 



t. Thus, we 



(16) 



Since we have not found a simple scaling form for global 
correlations in the kinetically constrained models, we test 
the hypothesis l|15[) directly. To this end, we chose pairs 
of total and waiting times t and t w such that the global 
correlation takes the same value, say C, 



C(tW,t£>)=C for k= 1,... 
and we plot the pdf's of local correlations 
P(C r (i«,4 fe >)) far * = !,..., 



(17) 



(18) 



The hypothesis l|15(l implies that for fixed value of C these 
pdf's should collapse onto a master curve. 

In Fig. \5\ we show the rather good demonstration of 
master collapse onto the Gumbel form for the modified 
FA model. To demonstrate that the strongly skewed dis- 
tributions that we observe arise from non-trivial dynam- 
ics induced by facilatation, we construct a toy model that 
contains the same lattice features as the FA model but 
evolves via unconstrained dynamics. This model is dis- 
cussed in detail in Appendix B. As shown in Fig. El the 
distribution of local correlations in this model is much 
weaker, and may be completely attributed to binomial- 
like fluctuations due to small coarse-grained cells. Fur- 
thermore, the data show no regime of approximate col- 
lapse onto a master curve, and for larger values of t w 
the distributions rapidly evolve to a Gaussian form. The 
non-trivial form and approximate scaling of the PDF of 
local correlations results from the dynamical constraints 
that induce a non-trivial dynamic length scale. 

Now we look at the d = 2, p = 2 version of the model. 
In Fig. {7\ we show the behavior of the normalized two- 
time connected correlation functions. As in the case of 
d = 3 p = 3, we could find no global scaling of the decay 
for all waiting times. Figured also shows the collapse of 
local correlations with C ~ 0.7 onto the Gumbel form. 
Interestingly, while the Gumbel parameters may depend 
sensitively on temperature, dimensionality and the value 
of C for which the collapse is plotted, the scaling onto the 
Gumbel form is remarkably accurate. 



C. Edwards- Anderson disordered spin system 

Similar results have been obtained for the disordered 
EA spin-glass model with usual Monte Carlo dynamics. 
In Fig.|Slwe show the pdf's of global correlations and in- 
tegrated responses for a small system. In this case, since 
the system is not macroscopic, we explore "mesoscopic" 
fluctuations that include similar information as the local 
fluctuations when L is not much larger than the corre- 
lation length. We see that the generic behavior of the 
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FIG. 5: Scaling of the normalized PDF of local correlations in 
the modified FA model, d = 3, L = 100, I = 5 and T = 0.4. 
C = 0.88 (top) and C = 0.70 (bottom). In both panels the 
data points correspond to five pairs of t and t w , with t w = 
10 , . . . , 10 time steps. The thin line is a Gaussian normal 
form. The dark solid lines are fits to a Gumbel form with 
parameter a = 11 (top) and a — 13 (bottom). 



data is very similar to the one displayed in the previ- 
ous subsection, although the particular time scaling [of 
the form I0J-I© in this case] is different. The data is 
also described for time differences that are shorter than 
the characteristic relaxation time by a negatively skewed 
form that resembles a modified Gumbel distribution. The 
trend in t — t w is the same one although the approach to 
the Gaussian form is much slower and this asymptotic 
form is hardly reached within the time window explored. 
We have also simulated the procedure used by Cipelletti 
et at, i.e. we have measured temporal fluctuations in 
mesoscopic samples, and found very similar results. 

It should be emphasized again that the glassy physics 
behind the EA model and the FA model differs in that 
a spin-glass transition exist at finite temperatures in the 
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FIG. 6: Scaling of the normalized PDF of local correlations 
calculated from the toy model for the parameters correspond- 
ing to the conditions described in Fig. Notice that the lack 
of data collapse in the dynamics toy model shows that the 
collapse we find in Fig. 5 is a consequence of a non-trivial 
correlated dynamics. 



EA model, while no thermodynamic transition exists at 
finite temperature in the FA model. The growing corre- 
lation length in the FA model is of a dynamical nature. 
We have verified that the Gumbel-like distributions sys- 
tematically revert to Gaussian distributions when the cell 
size is made larger. This suggests a direct way to estimate 
dynamical correlation lengths. Work along these lines is 
underway and will be reported in a future publication. 

VI. ANALYTIC ARGUMENTS 

In this section we recall and extend some of the the- 
oretical arguments introduced in22i2ii2£ to describe local 
and mesoscopic fluctuations in glassy systems. In those 
articles we focused on the time scaling of the PDF's with- 
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FIG. 7: Correlations in the modified FA model, d = 2 and 
L = 1000. Top: Scaling of the pdf's with C ~ 0.7 at the same 
temperature for mesoscopic blocks of size £ = 10. The data 
correspond to five pairs of t and t w , with t w = 10 , . . . , 10 5 
time steps. The dotted line is the Gaussian normal form. 
The dark solid line is a fit to a Gumbel form with a = 35. 
Bottom: Toy model results for conditions corresponding to 
the top panel. 

out paying special attention to the form of these distri- 
butions. Here, instead, we investigate their form and we 
develop an effective theory to describe them. 



A. Reparametrization invariance 

The two-time scaling of correlation functions that are 
monotonic functions of the two times take the form 

°^>~'(t§t) (19 » 

within each "correlation scale" as defined irJ°. In Sect. El 
we discussed the dynamic evolution of three kinetically 
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FIG. 8: Normalized pdf's of global correlation CV and inte- 
grated response \r (as defined irj22i) for a small 3d EA system 
with linear size L = 8, temperature T = 0.4 and waiting 
time t w = 1000. Data averaged over 8000 samples. The full 
line is the Gaussian distribution, while the dashed line is the 
generalised Gumbel form with a — n/2. 

constrained models and the 3d EA spin-glass. For the lat- 
ter, the aging decay is described quite well by a two-scale 
relaxation, see Eqs. J3J and (J5J, with a "simple aging" 
form, h$(t) = t, for the slow part. For the kinetically 
constrained models we have not found a single function 
h s (t) that describes the full slow decay from C = 1 to 
C = for all waiting times and times explored. 

In order to keep the presentation as simple as possible, 
we discuss the theoretical arguments that allow one to 
predict the behavior of local fluctuations, assuming that 
the full slow decay is represented by the scaling given in 
(|19fl with a single function k(t) — h s (t). By the end of 
this section we argue how this approach can be extended 
to describe the case of a more complicated slow relaxation 
that does not fall into this category. (The failure may be 
due to the development of many correlation scales in the 
slow relaxation.) 

In£^ it has been shown that the dynamic action con- 
trolling the local coarse-grained dynamics at long wait- 
ing time t w and long time differences t — t w in the EA 
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model approaches a form that is invariant under global 
reparametrizations of time that act as 

* -► h(t) , (20) 

with h(t) a monotonically growing function. This change 
transforms the "fields" in this action in a special way 
that, at the level of expectation values, corresponds 
to leaving invariant the local correlations, C' r (t,t w ) — > 
C r (h(t), h(t w )), and the integrated linear responses, 
Xr(t,t w ) —>Xr(h(t),h(t w )). Symmetry breaking terms 
exist but they become smaller and smaller (irrelevant) 
as times and time differences increase. These symmetry 
breaking terms (which have their origin in the details of 
the short-time and short-time-difference dynamics), al- 
though scaling down to zero at long times, are respon- 
sible for selecting a particular parametrization for which 
the global C decays from q EA to as t is increased. No- 
tice that if reparametrization invariance is left unbroken, 
C must be a constant independent of the times. Without 
entering the details of this action, let us introduce the 
main consequences of the existence of this approximate 
invariance. 

Due to the global invariance, and the scaling in (|19fl . 
one can argue that the slow part of the coarse-grained 
local correlations scale in time as22i2ii2£: 

c '<»»>~Kw)' (21) 

where the functions h r (t) are local time reparametriza- 
tions, t — > h r (t). The superscript S means that we are 
looking here at the slow part of the decay of the local 
correlation only (in other words, the decay for long time 
differences compared to the waiting time t w ). The argu- 
ment in i|21|) varies from unity when the times tend to be 
equal, t — ► t+, to zero at widely separated times t^$> t w . 
The external function / takes values between /(l) = q EA 
and f(0) = 0. The first limit is independent of the po- 
sition r since one does not expect the plateau to be a 
fluctuating quantity if a sufficiently large coarse-graining 
volume is used. The second limit holds if there is no 
external field applied and all local correlations decay to 
zero. We consider here monotonic local reparametriza- 
tions of the times, such that h r {t w ) < h r (t) at all r for 
t w < t. This is equivalent to assuming that the time lags 
we shall consider are not extremely short. 

Our original claim was that the external function / 
is the same for all coarse-grained centers in the sample, 
while all spatial fluctuations are encoded in the inter- 
nal function h r . The reason for this proposal is that 
the global reparametrization invariance in time of the dy- 
namic action in this two-time regime leads to low action 
excitations (Goldstone modes) for smoothly varying spa- 
tial fluctuations in the reparametrization of time, but 
not in the external form of the correlations. As in a 
sigma model, the external function / fixes the manifold 
of states, and the local time reparametrizations corre- 
spond to fluctuations restricted to this fixed manifold of 



states. The relation to the sigma model is described in 
detail in RefiSi and we shall not repeat it here. 

The symmetry-based argument and its predic- 
tions, originally developed on the basis of the time- 
reparametrization invariance of the dynamical action of 
the EA model, was proposed to hold beyond this model 
whenever a separation between fast and slow decay in the 
global correlation and response develops in time. Our re- 
sults on models without quenched disorder conclusively 
demonstrate the generality beyond the EA model. 



B. Effective random manifold action 

As in the sigma model approach to an interacting sys- 
tem in which we have identified a relevant variable and a 
relevant symmetry, we propose an effective action for the 
relevant degree of freedom that quantifies the fluctuations 
about the globally symmetric result. The philosophy is 
the same as the one followed when describing spin- wave 
excitations in a Heisenberg ferromagnet with a quadratic 
action that depends on the angular variables only. 

The scaling form in Eq. (|21|l can be written in an equiv- 
alent form by defining 4> r (t) = hih r (t): 

C s r (t,t w ) ~ /(e-t-M*)"^)]) 

= /(e- A ^lt.) . (22) 

In the dynamic case under study, after parametrizing the 
correlation in this way, the relevant field is the function 
4> r {t). The action governing the dynamics of 4> r {i) de- 
pends on the details of the particular problem. However, 
we can greatly restrict the form of the possible actions 
by simply using the constraints due to the symmetries. 
These are: 

i. The action must be invariant under a global time 
reparametrization t — > s(t). 

ii. If our interest is in short-ranged problems, the ac- 
tion must be written using local terms. The action 
can thus contain products evaluated at a single time 
and point in space of terms such as 4> r (t), <fi r {t), 
V0 r (t), V</> r (t), and similar derivatives. 

Hi. The scaling form in Eq. 12 I'll is invariant under 
4> r {t) — > (j> r {t) + with <i> r independent of time. 
Thus, the action must also contain this symmetry. 

iv. The action must be positive definite. 

These requirements largely restrict the possible ac- 
tions. The action with the smallest number of spatial 
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derivatives (most relevant terms) is 



d d r / dt 



K 



(v<M*))' 

4> r {t) 



M (j) r (t) 



(23) 



The last term is simply a total derivative, and its space- 
time integral is constant, so we shall drop it. Hence, we 
find that the strong constraints imposed on the action 
reduces it to the form 



SV.ff =K d d r dt 



4>r{t) 



(24) 



Notice that the action now solely depends on the time 
derivatives <f> r (t). 

Let us now examine the consequences of this action 
on the form of the local correlations Ij22(l . Due to the 
simple form l|24(l the <p r are uncorrelated at any two dif- 
ferent times t\ and ti- Thus the expression A0 r |( m = 
f t dt'(j) r (t') entering the exponential in the scaling form 
in Eq. 1221) is a sum of uncorrelated random variables 
in time. Hence, one can interpret such expression as 
the displacement of a random walker with position de- 
pendent velocities. Alternatively, one can think of the 



space-dependent differences A(j) r 



/. dt'(j> r (t') as the 



net space-dependent height (labeled by t) of a stack of 
spatially fluctuating layers dt'<fi r (t'). The action for the 
fluctuating surfaces of each layer is given by Eq. I|24(l . 

Before proceeding, let us take a moment to show 
that indeed the action in Eq. (|24|l satisfies all the four 
constraints enumerated above. First, consider a global 
change of variables t — > s(t), for which the action trans- 
forms as 



K d d r dt 



K d d r ds 



(f ) 2 (vjrKg) 

m Mm) 

(v<Ms))' 



s E 



(25) 



Hence, as we required in point i., the action is invariant 
under the global time reparametrization. The action is 
also clearly local. Because it only contains explicitly 4> r 
but not 4> r , it is invariant under 4> r {t) — > (j) r (t) + $ r , 
with <£v independent of time. Finally, since the lo- 
cal reparametrizations h r (t) are monotonically increas- 
ing functions of time, so are <j) r {t) — lnh r (t), and conse- 
quently 4> r (t) > 0. Hence, the combination (V 4> r ) 2 / 4> r is 
always positive, and the action is positive definite. 

To better understand the consequences of this effec- 
tive action, we study it in two regimes. The first is one 
in which we consider only small spatial fluctuations of 
the local reparametrizations with low velocities. In the 
second, we consider the full effect of the action, without 



assuming that the local fluctuations are small. We will 
use the latter to obtain the full form of the probability 
distribution of local correlations and compare it to the 
numerical and experimental results. 



1. Linear regime 

We begin by making some simplifying assumptions so 
as to obtain a better understanding of the consequences 
of local reparametrizations22i2£. We shall shortly relax 
these assumptions. If we assume that the site-to-site 
fluctuations are small, the field <fi r (t) will be given by 



(f> r (t) = s(t) + ip r (t) with ip r (t) < s(t) , (26) 

and the uniform part s(t) defined as 

s(t) = ]nh s (t) (27) 

with h s (t) the scaling function of the global correlation. 
In terms of the s parametrization we can write 



</> r (s) = S + <f r (s) , 



and 



(28) 



(29) 



4> r (s) = 1 + tp r (s) . 
We shall further assume that 

(fr{s) < 1 . (30) 

This approximation leads to the action 

S EFF J d d r J ds {Vifr(s)f • (31) 

The surfaces tp r (s) are uncorrelated for different 
reparametrized times s. The argument that enters in 
the evaluation of the local correlation functions C r (t, t w ) 
is, according to Eq. I|22l) . 



(O 



ds [1 + <p r (s)] 



= s(t) - s(t w ) + \J s(t) - s(t w ) X r 



h s {t w ) 



h s (t w ) 



(32) 



where the square root arises from the fact that the <p sur- 
faces are uncorrelated for different s, leading to a sum of 
independent variables in the reparametrized time direc- 
tion. The X r is a Gaussian random surface that encodes 
the spatial fluctuations, and it is distributed with prob- 
ability 



P(X r )oce~ K I ddr {vx ^ 2 



(33) 
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Even within the linear approximation l|29|) and <|30|) . 

one can already explain why the probability distributions 
for local correlations of a system with a global scaling as 
in iQJ-© collapse as a function of h s (t) / h s (t w ) . This fact 
is a consequence of Eqs. l{2"2"|) and (J23! 20 '3i. Moreover, 
noticing that the two-time global correlation is in one- 
to-one correspondence with the ratio h s (t) / h s (t w ) , one 
concludes that the PDF of local correlations scales in time 
as the global correlation itself. 

While this approximation gives the observed time scal- 
ing of the PDf's of local correlators, it does not give their 
correct shape, in particular when t — > <+. The reason 
is that the assumption that site-to-site fluctuations and 
their velocities are small as compared to the bulk average 
fails more severely in this regime. As we show next, we 
can still tackle the problem in the non-linear regime. 

2. N on- approximate treatment of fluctuations 

In order to get a more accurate prediction of the form 
of the PDF, one needs to go beyond the approximations 
in Eqs. (|26[1 and J2UJ. As done in the previous section, 
we start by writing the action in the parametrization 
s(t) = In h s (t) chosen by the bulk values, and by defining 
a new fields, tpr(s), such that ip1{s) = 4> r (s), one then 
has 

S EFF = K j d d r J ds (Wv(s)) 2 . (34) 

This is again the action of uncorrelated Gaussian surfaces 
for different proper times s. 

The argument entering in Eq. (1221 now reads 

,s(t) 

A0 r |* iu = / dsMs)= / ds<ip*(s), (35) 

Js(t„) J s(t„) 

Due to the Gaussian statistics of the Vv(s), it is simple 
to show that connected TV-point correlations of A<j) ri \\ m 
satisfy 

(A0 r X A<f> r2 \l ■■■ A0 r X> C = ( 36 ) 
[s(i) - s(t w )] T(r ll r 2l ■ ■ - ,r N ) , 

where the function T can be obtained from Wick's theo- 
rem, summing over all graphs that visit all sites (con- 
nected) with two lines (because of ip 2 ) for each ver- 
tex i corresponding to a position r^. Notice that the 
reparametrized times appear only in the prefactor As = 
s(t) — s(t w ). Therefore time dependencies are functions 
of As alone. Once again we find that the probabilistic 
features of the fluctuations of local correlations are func- 
tions of the scaling variable, As — In h s (t) / h s (t w ) , and 
hence of the global correlation itself. This coincides with 
the numerical observations. 

Here, a comment is in order. In obtaining the effec- 
tive action, we have neglected higher gradients and time 
derivatives because these were considered irrelevant. In 



particular, this is justified for large reparametrized time 
differences. Note, however, that these terms are im- 
portant when As becomes of order unity. This short 
(reparametrized) time-cutoff scale means that the sur- 
faces for different si,2 are uncorrelated only when |si — S2 1 
is larger than a cutoff Ss. Keeping this in mind, let 

n 

Acp r \l = Ss V£(*0 . ( 37 ) 

i=l 

where we have discretized the integral in Eq. I|37|) in steps 
of the cutoff 5s. The number of such steps is n — As/<5s. 
The expression (|35ll and the scaling relation l|22[) allow 
us to compute the PDF for the local correlation functions 
(these expressions can also be used to investigate TV-point 
correlations of C r (t,t w )). We shall do this in the next 
section generating numerically random Gaussian surfaces 
■i/v(s), from which we obtain A</v|^ and then C r (t,t w ). 
Let us now discuss the extreme cases (when n is either 
small or large) before turning to the numerical calcula- 
tions, which are better suited to describe the crossover 
between these two limits. 

Case C ~ q EA 

If the global correlation is close to the Edwards- 
Anderson parameter, n — As/Ss is small. Consider for 
simplicity n = 1. In this case, letting Y r = ip r (si) 1 we 
write 

A0 r \l =5sY r 2 , (38) 
with Y r distributed according to 

P(Y r ) OC C- K S d d r _ (3g) 

Using Eq. (1221) the local correlation is given by 

C r (t,t w ) = f{e- A ^L) =/(e-^ 2 ) . (40) 

In what follows, we connect the PDf's of the local corre- 
lations C r coarse grained over a volume around the point 
r and the PDf's for the fluctuations of the global mag- 
netization m of the xy magnet at low temperature, as 
studied by Bramwell et alm^, or the roughness of a Gaus- 
sian surface of the Edwards- Wilkinson type, as analyzed 
by Racz et alml. For concreteness, we use in the following 
the language of the XY model. 

In the case of the xy magnet, the local field 8 r is dis- 
tributed according to 

P(6 r ) cx e~ K f ddr ( ve "> 2 , (41) 

and one can write the local magnetization in one direction 
as 

m r — cos 8 r , (42) 

where for simplicity we took m = 1. Bramwell et al. 
found that the distribution of the global magnetization 
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of finite-size systems does approach a scale- invariant non- 
Gaussian distribution (when the magnetization is nor- 
malized by subtracting the average and dividing by the 
standard deviation). In particular, the PDF is skewed and 
can be approximated by a generalized Gumbel distribu- 
tion. The underlying reason for the non-applicability of 
the central limit theorem is that the system is critical 
(there is no finite correlation length) and hence the local 
magnetizations are not uncorrelated. Similarly in a sys- 
tem where the correlation length £ is finite, non-Gaussian 
distributions will be found as long as the coarse graining 
length is smaller than £. 

The connection between the PDF's P{C r ) and P(m) 
for dynamic local coarse grained correlations C r (t,t w ) in 
our case and the global magnetization m in the xy model 
is as follows. First, the distribution of the variables Y r 
and 6 r in Eqs. i|39|) and (|41|l are identical. The fact that 
the functional dependencies of C r and m r on Y r and 9 r 
respectively, as given in Eqs. 14U|) and H42|) . are not the 
same lead to slightly different forms of the PDF's that, in 
both cases, can be approximated by generalized Gumbel 
forms with different parameters. To illustrate this point, 
consider the case in which the stiffness K is large so that 
both variables Y r and 9 r (centered around zero) are small 
within a given coarse graining volume. In this case, we 
can expand Eqs. (|4L)I) and 14211 as follows: 



C r = f(e- 5sY ^ Y r 2 /'(1) + .. . (43) 

m r = cos 6 r w 1 - - 6l + . . . (44) 



and because the PDF's are calculated for normalized vari- 
ables (zero average and unit variance), the constants in 
the problem (such as /(l) and /g(l)) have no effect on 
the distributions. Hence, in this case the PDF's for the 
two problems are identical. The sum over a finite vol- 
ume, used to compute the global mesoscopic m does not 
modify this conclusion. In particular, in 2c? the PDF can 
be approximated by a generalized Gumbel distribution 
with parameter a « i/2 (see App. 0, as discussed in 
detail in Refi2£. 

We want to stress that the description of the PDF's 
with Gumbel-like extreme statistics is, in general, just 
an approximation that correctly captures their main fea- 
tures, namely a single maximum, negative skewness that 
crosses over to a positive skewness, and an exponential 
or near exponential tail for small fluctuations below the 
mean. Moreover, C r is not a priori an extremal quantity 
so there is no deep reason why it should be distributed 
following the extreme value statistics. 

Case C g 

If the global correlation is close to zero, n = As /5s is 
large. In this case, the sum in Eq. (|37() has many terms, 
and since all terms in the sum are positive, the total will 
have some average that is proportional to As plus fluc- 
tuations. Some reflection point us to the same situation 



that we presented in Sect. I VI SH where we discussed the 
linear regime. Thus, we can write the local correlations 



/As Xr 



(45) 




where the variable X r is distributed according to 
Eq. i|33|) . From here one can show that if f(x) is a mono- 
tonic increasing function, the PDF P(C r ) is positively 
skewed (as opposed to what is found when C rs q EA ). 
Indeed, expanding C r around X r ~ 0, 



Cr 



one sees that the signs of the quadratic term in X r in 
Eq. i|T7|) and the quadratic term in Y r in Eq. have 
opposite signs. Since this sign dictates the skewness (re- 
call that X r has a symmetric distribution) , one concludes 
that the PDF is positively skewed in this limit, as also ob- 
served numerically and experimentally. 



C. Beyond a single scale 

The argument used above can be heuristically ex- 
tended to deal with more general problems. Let us as- 
sume that global correlations evolve in a sequence of 
scales, each characterized by their own scaling function, 



hk(t), and external function, /; 



c(t,t w ) = j2fk 
fc=i 



23,50.57 



hk(t w ) 



h k (t) 



(47) 



with each hk (f ) different monotonically growing function 
of the original time. This means that when the ratio for 
a chosen value of k varies as, 



o < hk(tw) < 1 , 



h k {t) 



all others are constant: 

hk>(<k){tw) _ q 



hk'(>k){tw) 



1 



hk'(<k)(t) fyfe'(>k)(*) 

The boundary conditions are such that 

lim f k (x) = q k , lim f k (x) = 



(48) 



(49) 



(50) 



with < qk+i < qk < 1, and J2k =1 qk = 1- In the time 
regime in which the k-th ratio varies the global correla- 
tion reads 



C(i, t v 



fk 



hk(t w ) 
h k {t) 



«*' 

k'>k 



(51) 
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Motivated by what was presented mSLSii we propose that 
the local fluctuations can be parametrized as 

<««.<■> -tA(W 

n 

= J2h(^ kAtm) ' 4 ' kAt) ) ■ (52) 
k—1 

We assume that the local fluctuations are small in the 
sense that they do not exchange scales. This is consistent 
with the supposition that the Edwards- Anderson param- 
eter, or more generally in this case the plateau values q^, 
are not modified by the fluctuations. 

Based on the symmetry arguments already explained, 
we propose that each (j>k r (t) is distributed as in This 
leads, after the reparametrization of time, to 

s k (t) = hih k (t) , (53) 

and the change of variables (ipk r ( s k)) 2 — 4>k r (t)> to a 
Gaussian distribution for each ipk r {sk), 

St F =K j ' d d r J ds k (VW(sfc)) 2 ■ (54) 

Note that the </>fc r (sfc) are positive and we assumed that 
the stiffness K is k independent, for simplicity. This is 
simply due to the fact that the hk (t) characterizing each 
two-times sector are monotonically growing functions of 
time. 

Thus, when we look at the connected correlations be- 
tween the A(j)k r \tu, ' f° r given values of t and t w that corre- 
spond to the scale k for the global correlation, only the k- 
th scale makes a non-trivial (non-constant)contribution. 
Since in each two-time sector the global scaling of the 
PDF of local correlations is the same as the one for the 
global C, wc find that this result holds for all times. Once 
we have analyzed the case of n scales, we can take the 
limit n —> oo that corresponds to the ultrametric scaling 
form and claim that in this case the pdf's scale as the 
global ones. Thus, we expect the collapse of the type 
found for the facilitated models. 



D. Numerics for random surfaces 

In this section we numerically generate random sur- 
faces with the statistics given by (|35|l and we calculate 
the pdf's of local correlations related to the surfaces us- 
ing Eq. (|22[1 . In order to carry out numerical evaluations 
of the local C r (t, t w ) pdf's, let 

n 

Act> r \\ w = 5s , (55) 

i=l 

where we have discretized the integral in Eq. i|35|) in steps 
of the cutoff 8s. We proceed as follows: 



i. We generate independent Gaussian random sur- 
faces ip T {si), which we choose to be of unit variance 
since the choice of variance does not change the re- 
sults qualitatively. These surfaces are generated 
with the distribution P[tp r (si)} oc e -K ! d A r{vi> T f 
(see Eq. ®). 

ii. We then sum the ^r(sj) over n time slices. The 
number of slices determines the time difference t — 
t w . 

Hi. Next we generate the surfaces C r (i, t w ) — e~ A ^ r ^^ , 
i. e. we choose for simplicity the "external" function 
f(x) = x [see Eq. (QJ]. 

iv. Finally, we look at the pdf's P(C r ) of the nor- 
malized coarse-grained local correlations C r (sub- 
tracting the average, scaling by the standard de- 
viation, and coarse-graining over a box of volume 
i 3 = (2M + l) 3 ) as a function of the number of 
slices n. 

Before presenting the numerical results, let us discuss 
in more detail the issue of coarse-graining the correlation 
C r over a certain volume. There are actually two levels 
of coarse-graining for theories defined on a lattice. First, 
in order to define continuous fields h r (t), cf> r (t), and ulti- 
mately ip r (t), a coarse- graining of a lattice theory has al- 
ready occurred. Coarse-graining the local correlation C r 
over a given volume is a second level of coarse-graining. 
The connection to the XY model and the fluctuations of 
the global magnetization studied by Bramwell et al£& is 
helpful in clarifying this point. Starting from a lattice 
theory, one coarse-grains to obtain a Landau-type the- 
ory of a continuous variable m r . Then, to study the fluc- 
tuations of the total magnetization of the system over a 
given volume, one must average m r over a second coarse- 
graining volume (which in the case of the global magneti- 
zation is the entire volume of the system) . In the numer- 
ical simulations we present here, we start from an action 
for the tp r that is Gaussian. The first level of coarse- 
graining that led to an action for a continuous 4> r or ip r 
was assumed to have taken place already. So even though 
in the numerical calculations we employ a discretization, 
one should keep in mind that a first coarse-graining has 
been invoked. 

In Fig. El we show for a system of size L = 32 the de- 
pendence of the P(C r ) on the number of slices. We show 
curves for n — 4, 8, 12, 16, 20, and for a coarse-grained 
volume £ — 3. The smaller the value of n, the smaller 
is the ratio h s (t)/h s (t w ). Several comments are in or- 
der. First, notice that for small n, the curves clearly 
deviate from a Gaussian form and are reasonably well 
fit by a generalized Gumbel distribution (see App. A). 
The curves become more Gaussian as n increases. This 
behavior is qualitatively similar to the one we obtained 
for the facilitated spin model and the 3d EA model we 
studied numerically in Sect. as well as the experimen- 
tal observations of Cipelletti et altMl. The approach to a 
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(C r -<C r >)/a 




-6 -4 -2 2 4 



(C r -<C r >)/G 

FIG. 9: Normalized pdf's P(C r ) for local coarse-grained cor- 
relations C r generated by the random surfaces suggested by 
the analytical theory. The curves shown are for sums over n 
slices. As n increases they become less skewed and approach 
the Gaussian distribution (full line). Top: d — 3, system of 
linear size L = 32, and the coarse-graining cell of linear size 
£ = 3. The Gumbel parameter is a = 2.2, 4.5, 11 (for increas- 
ing n). Bottom: d = 2, L = 32, and £ = 5; a = 0.4, 2, 10 (for 
increasing n). 

Gaussian behavior is however very slow and we do not 
reach the Gaussian form with the values of n used. 



E. Essential ingredients for universal fluctuations 

Why do the analytical arguments we presented above, 
based on a symmetry (reparametrization invariance), 
lead to the universal behavior and collapse (for a fixed 
global correlation) of the P(C r ) that is observed in the 
facilitated spin models, in the 3d Edwards-Anderson 
model, and in some experiments^? We can single out 
two essential elements of the theoretical description that 
allow it to account for this behavior. 

The first ingredient is that the theory proposed in 
j^ e f S| 20 i 2i i 3^ contains a massless Goldstone mode as- 



sociated with the spatial fluctuations in the asymp- 
totic long-time (and long-time-difference) limit. This 
means that correlation functions such as those in 
Eq. (|37|) are described by algebraically decaying functions 
.F(ri, T2, ■ ■ ■ , rjv)- The absence of a correlation length 
£ leads to the non-Gaussian behavior for quantities av- 
eraged over a finite volume, similarly to the studies in 
Refsi2£i2£. Of course, one expects that universal behavior 
of the pdf's will occur only after coarse-graining. Note 
however that it is the absence of a correlation scale that 
implies that the law of large numbers will not take ef- 
fect as the box size increases. So the theory based on 
reparametrization invariance of the action gives rise to a 
non-Gaussian distribution for coarse-grained C r (t, t w ) (in 
addition to any coarse-graining needed to define a theory 
in the space continuum) because of the algebraic correla- 
tions originating from the massless Goldstone mode. In 
the theory of Refsi^2*2ii2£ we argued that the correlation 
scale £ — ► oo in the asymptotic long-time limit. However, 
for finite times there is a finite Coarse-graining much 
beyond this length scale leads to Gaussian distributions. 

The second element that the analytical theory contains 
is that all dependencies on the times are functions of 
As = s(t) — s(t w ). This leads to the collapse of distri- 
butions (not only their bulk averages) as a function of 
In h s (t) / h s (t w ) . In addition, the fact that the different 
n = As /5s slice contributions ip^.{si) are uncorrelated, 
and that they enter in the local correlations C r (t,t w ) 
through the sum over all slices, is the reason why the 
skewness of the distributions depend on n, as shown in 
Fig. El Again, this is in good accord with the experimen- 
tal observations of RefpiS and with the numerical results 
for the models that we study in this paper. 

In summary, we believe that the qualitative features 
shared by the pdf's of local correlators in the facilitated 
models, in the 3d EA model, and in the experimental 
observations of Cipelletti et alMi can be understood as a 
consequence of certain spatial and temporal correlations 
of the C r (t, t w ), which are contained in the theory based 
on reparametrization invariance. 



VII. CONCLUSIONS 

In equilibrium and away from criticality any global ob- 
servable of a macroscopic system has Gaussian fluctua- 
tions. At criticality, instead, one observes non-Gaussian 
fluctuations due to the divergence of the correlation 
length and the non applicability of the central limit the- 
orem. Still, scale invariance at the critical point con- 
strains the possible probability distributions; these are 
determined by the universality class to which the systems 
belong. A similar criterion to classify the probability dis- 
tributions of the fluctuations of macroscopic observables 
in critical noncquilibrium systems is based on the use of 
symmetries. It has been proposed in2£. 

The glass transition, where the system falls out of equi- 
librium, is a dynamic crossover and neither a dynamic 
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nor a thermodynamic transition. However, many fea- 
tures of glassy dynamics resemble critical relaxation such 
as the fact that correlations do not decay exponentially, 
but with much slower decay forms. One could expect 
then that some concepts that have been useful to study 
critical phenomena could also apply to the glassy dynam- 

In this paper we continued the study of local fluctu- 
ations in the dynamics of glassy systems. On the the- 
oretical side we improved the sigma model, or random 
manifold action, proposed in22i2iiS to capture not only 
the scaling in time of the PDF's but also their functional 
form. On the numerical side, we showed that part of 
the predictions of the theory tested inSLSl using the 3d 
Edwards- Anderson model also hold for a non-disordered 
kinetically constrained spin model of a glass. This con- 
firms that the existence of quenched disorder is not im- 
portant in this respect and that kinetic frustration and 
energetic frustration lead to very similar non-equilibrium 
dynamics. Moreover, we tested the evolution of the form 
of the PDF's for these two models and verified that they 
behave in a way that resembles strongly the experimental 
observations iniS. 

One could expect that a more detailed and extended 
analysis of the local dynamics of glassy systems might 
lead to better differentiation between models and sys- 
tems. For instance, one could find that the simple action 
(|35|l is not enough to reproduce the special form of PDF's 
found in a system and hence be forced to include other 
terms that we here neglected. 
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APPENDIX A: GENERALIZED GUMBEL 
DISTRIBUTIONS 

The Gumbel distribution is defined as 

* B (y) = ^(«<™>-«<-» >) , (Al) 

r(o) 

where T(a) is the gamma function. The parameters yo 
and a control the position of the center and the width of 
the distribution, respectively. One could define standard 
distributions by requiring that the center is at zero, and 
the width is unity. These conditions fix, for a given a, 
yo(a) and a(a) as 



a = y/W(aj , ay = logo- *(a) , (A2) 

where ^(a) is the digamma function ^(a) = T / (a)/T(a). 

The Gumbel distribution with a — 1 appears as 
one kind of asymptotic extreme value statistics when 
searching for the distribution of the maximum (or the 
minimum) of a sequence of independent identically dis- 
tributed random variables with a probability density de- 
caying faster than any power law. Extensions with inte- 
ger parameter a > 1 are the results of searching for the 
distribution of the a-th largest (smallest) value in the se- 
quence. For any a, there is a choice of sign for a and the 
two signs correspond to the Gumbel statistics of cither 
extreme minima or extreme maxima. Further extensions 
with non-integer parameter a < 1 have been found when 
the elements of the sequence are correlated over a dis- 
tance l/a£&. For a = tt/2 one recovers the distribution 
studied by Bramwell, Holdsworth and Pinton 59 . 

Finally, let us study two limits of the generalized 
Gumbel-like distribution (|A1|) . A Gaussian distribution 
of unit variance is obtained by taking the a — > oo limit 
of $ a (y) while holding the product a 2 a — 1 fixed: 



lim $ a (y) 

a — *oci.a A a— 1 



lim 



a^oo,a 2 a=l y/a T(a) 



a a(a(y-y a )-l-a(y-y )-±a 2 (y-y ) 2 +O(a 3 )) 



lim c -Uv-vo) 2 

a^oo y/a r(a) 

1 



e -^(y-yo) 2 



(A3) 



where in the last step we used Stirling's formula to get Another interesting limit is when a — > while holding 
the normalization factor. \a\a = 1: 
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lim <!> a (y) = I™ e sgnQ 9 (-sgna (y - y )) 

D,|a|a=l a-»o al (a) 

= c sgna 9 (-sgna (y - t/ )) (A4) 



where we used that lim a ^ aT(a) — T(l) = 1. Thus, 
the a — > limit leads to exponential distributions either 
starting or ending at y = yo. 

APPENDIX B: INTRINSIC SKEWNESS OF 
UNCORRELATED DYNAMICS 

The probability distribution of local correlation func- 
tions evaluated at finite times and coarse-grained over a 
sufficiently large cells is, to a first approximation, nor- 
mally distributed around the bulk value. However, for 
the times and cell sizes considered in this paper, one 
needs to consider how the intrinsic distribution of uncon- 
strained spins contributes to the skewing of the distribu- 

I 



Here n is the number of spins up initially up, c is the ini- 
tial bulk concentration of up spins, d the number of ini- 
tially up spins that are down at time i, and u of initially 
down spins which are up at time t. In the equilibrium 
case, p and q are related by 

Qeq = c eq /(l - c eq )p eq (B2) 

The same relationship (without the subscripts) holds for 
a constant concentration system, but not in the case of 
a steadily decreasing concentration, as in the quenched 
facilitated models. That however can be resolved by in- 
cluding a scaling factor X proportional to the concentra- 
tion during the quench as a function of t w and t. 

For a given ensemble average value C of the single site 
correlation, one can obtain the PDF for the correlation 
averaged over a number V r — £ d of uncorrelated spins. 
(Notice that, since these are free spins without any dy- 
namical constraint, the number of time steps necessary 
to reach the average value C is much smaller than that 
for the kinetically constrained models.) Using the result 
in (|Blfl . one can obtain the corresponding value of the 
cell connected correlation function 

C c , norm (t) = c(0)c(i) - c^O) x tffj (B3) 

where c(i) = c(Q) = § , and c(0)c(t) = The 

overline refers to the sub-ensemble average. 



I 

tion of local correlations. In this Appendix, we consider 
the degree of intrinsic skewing that arises from uncorre- 
lated dynamics on a lattice. 

To differentiate between correlated and intrinsic skew- 
ness, we look at a model of paramagnetism without any 
facilitation. The Hamiltonian is then the same as in the 
FA model. This toy model has no dependence on the 
dimensionality or the neighboring environment, and de- 
pends only on the number of spins included in the cells. 
Therefore, all sites are independent and the probability p 
that an initially up spin is down at time t is the same ev- 
erywhere {q is the equivalent probability for the spins ini- 
tially down). Using Metropolis Monte Carlo single spin 
flips, one finds the following distribution for a cell of N 
spins using binomial distribution arguments, 



(Bl) 

I 





10 1 


10* 


10 a 


10 4 


10 b 


t 


30 


45 


90 


330 


1000 


c 


0.30865 


0.27257 


0.21203 


0.15505 


0.1685 


X 


0.938 


0.969 


0.988 


0.995 


0.9995 


P 


0.233 


0.230 


0.237 


0.258 


0.264 


a 


100 


~ oo 


~ oo 


~ oo 


~ oo 



TABLE I: Unconstrained parameters for a decay of the global 
correlations to a value of 0.7 of the initial value for the FA-2d2n 
model for a quenching T = 0.6 and a box size of AT = 10 x 10. 
The parameters are the t w defect concentration c, the concen- 
tration scaling factor X for time t, and the probability p of 
defect disappearance at time t. The resulting approximative 
Gumbel a parameter is also given. 



Using bulk values for the different parameters, the sta- 
tistical component of the distribution of the local spin 
correlation can be obtained, by numerically evaluating 
Eq. (|B3(1 for all n, u, and d and properly weighting them 
according to Eq. I|B1|I . Tables HI ITT1 and IIIII show the pa- 
rameters used in the different models. In the case of the 
plaquette models, the same discussion applies, but with 
a null magnetic field instead. 

In the first two cases (Table [I] and [HJ the bulk correla- 
tion is chosen to have a rather large value and the Gum- 



V(n, d, u) 



Nl 



{n - d)\d\{N - n - u)\u\ 



l (l - c) N - n p d (l - p) n - d q u (l - q) 



N-r 



19 



till 


10 1 


10* 


10 3 


10 4 


10 b 


t 


15 


26 


81 


370 


2300 


c 


0.3299 


0.2903 


0.2399 


0.2007 


0.1685 


X 


0.959 


0.982 


0.993 


0.996 


0.998 


P 


0.220 


0.219 


0.2305 


0.241 


0.249 


a 


80 


80 


100 


~ oo 


~ oo 



TABLE II: Unconstrained parameters for a decay of the global 
correlations to a value of 0.7 of the initial value for the FA-3d3n 
model for a quenching T = 0.8 and a box size of N = 5 x 5 x 5. 
Parameter definitions are identical to those in Table U] 



bel parameters found are extremely large corresponding 
to approximately Gaussian distributions. No interesting 
effect is observed in these cases. In the last example con- 
sidered (Table IIII|) the bulk correlation is chosen to take 
a smaller value and the parameters found indicate that 
the distribution is skewed. However, if one compares to 
the results obtained for the same model with facilitation 
one concludes that the non-Gaussian effect is here much 
weaker and tends to disappear rather quickly when the 
times involved get long. 



til] 


10 1 


10* 


10 3 


10 4 


10 b 


t 


2 


4 


10 


32 


100 


c 


0.3299 


0.2903 


0.2399 


0.2007 


0.1685 


X 


0.992 


0.997 


0.9996 


1.000 


1.000 


P 


0.0812 


0.0908 


0.0902 


0.0962 


0.0993 


a 


15 


20 


40 


90 


100 



TABLE III: Unconstrained parameters for a decay of the 
global correlations to a value of 0.88 of the initial value for 
the FA-3d3n model for a quenching T = 0.8 and a box size of 
iV = 5 x 5 x 5. Parameter definitions are identical to those in 
Table |H 



In summary, the simple correlations that follow from 
coarse-graining of some binomially distributed single-site 
spin-spin correlation whose average is restricted to be C 
cannot account for the PDFs we find for the kinetically 
constrained model. Spatial correlations between the spin 
variables are needed to account for both the width and 
skeweness of the distributions. 
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